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We reformulate Fourier-space crystallography in the language of cohomology of groups. 
Once the problem is understood as a classification of linear functions on the lattice, restricted 
by a particular group relation, and identified by gauge transformation, the cohomological de- 
scription becomes natural. We review Fourier-space crystallography and group cohomology, 
quote the fact that cohomology is dual to homology, and exhibit several results, previously 
established for special cases or by intricate calculation, that fall immediately out of the for- 
malism. In particular, we prove that two phase functions are gauge equivalent if and only if 
they agree on all their gauge-invariant integral linear combinations and show how to find all 
these linear combinations systematically. 
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The discovery in 1984 of crystals with five- fold symmetry, and therefore no periodic- 
ity, revived interest in Bienenstock and Ewald's 1962 reformulation of crystallography in 
Fourier space. ^ Mermin and collaborators have applied this Fourier-space crystallography 
to the classification of the space groups of quasi-periodic and periodic crystals, modu- 
lated crystals,^ and to color groups."^ We show here that this technique, by now familiar 
to crystallographers, has a particularly simple expression in the language of homological 
algebra. 

Section 1 quickly reviews the basis of Fourier-space crystallography. The phase function 
$g(k), defined for k in the reciprocal lattice and g in the point group, facilitates the classifi- 
cation of space groups. We review the gauge equivalence of two phase functions. In section 
2, we consider the necessary extinctions that occur in the diffraction pattern when g leaves k 
invariant, yet ^^(k) does not vanish. Life would be simple if two phase functions that agreed 
at all such pairs {g, k) were gauge-equivalent, but in fact there are two crystallographic coun- 
terexamples; however, when we consider gauge equivalence of linear combinations of phase 
functions, we can state the result (Theorem 1) that two phase functions are equivalent if 
and only if all their gauge-equivalent linear combinations agree. Section 3 reviews work by 
Mermin and Konig on the most important such linear combination, one that relates to a ray 
representation of the point group and necessitates an electronic degeneracy in the system. 
Section 4 places the phase function $ in a cohomology group, while section 5 puts gauge- 
invariant linear combinations of phase functions in a related homology group. In section 6, 
we invoke the duality of these cohomology and homology groups to prove Theorem 1 (an 
alternative, elementary, proof being presented in the appendix). We also prove that for any 
phase function, there exists a gauge in which it takes only rational values (with set denom- 
inator). Section 7 applies our results to simplifying the classification of space groups, and 
Section 8 addresses the ray representation of Section 3 in homological language. Section 9 
compares our results to a real-space treatment. 

1. Fourier-Space Crystallography 

Mermin and others have argued persuasively^'^ that the theoretical significance of qua- 
sicrystals lies not so much in relaxing the requirement of periodicity as in replacing exact 
identity of a density function {e.g., electronic or nuclear) under symmetry operations with 
indistinguishahility. We will review Fourier-space crystallography only tersely, in order to 
establish the notation; for a more developed exposition, we refer to Ref. [10]. 

The reciprocal lattice, C, (or simply lattice, since there is no direct-space lattice 
of translations in the aperiodic cases) is the closure under addition of the set of Bragg 
points k at which the Fourier transform p(k) of the direct-space density function does not 
vanish. Crystals, periodic or quasiperiodic, have reciprocal lattices generated as all integral 
linear combinations of a finite set of vectors k; we do not consider non-crystals. A different 
density function, p' , is said to be indistinguishable from p if all spatially- averaged n-point 
correlation functions, representing macroscopic physical measurements, are the same for the 
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two densities: in Fourier space, 

p'(ki)p'(k2) . . . p'{K) = P(ki)p(k2) . . . p{K) (1) 

for all lattice vectors summing to zero (J^kj = 0). This holds if and only if 

p'(^i,) = e^-ixil^)p{k) , (2) 

where called a gauge function, is any real-valued linear function (taking values modulo 
1) on C We are particularly interested in the case that the first density function in (2) 
results from acting on the second with an element 51 of a point group, G C SO (3), for then 

p'(k) = p(^k) = e2-^^(k)^(k) , (3) 

where the particular gauge function for this g, is called a phase function. Phase 
functions are constrained by the group-compatibility condition for g,h & G, 

^gh{k)^^g{hk)+^h{^) , (4) 

a consequence of the associativity of the group action on the lattice. 

The phase function defines algebraically the action of the point group on the density and 
so, together with the lattice and orientation, encodes the information of the traditionally- 
defined space group. However, altogether too much freedom remains, for two phase func- 
tions, and related by a gauge function through 

$;(k)-$,(k)=x(^k-k) (5) 

yield indistinguishable densities. To classify space groups, one therefore first computes all 
possible phase functions satisfying (4), then groups into gauge-equivalence classes phase 
functions differing only by a gauge, as in (5). 12 Until now, the second step has required a 
clever choice of specific gauge in which all but a few values of $g(k) vanish; homological 
algebra provides a more systematic approach to the classification, for it it is precisely the 
group of functions satisfying (4) with the identification (5) that constitutes the fundamental 
object of cohomology. 

2. Extinctions, an untrue proposition, and a theorem 

The class containing the trivial for which all values $g(k) vanish, is called symmor- 
phic. Perhaps the most direct result of Fourier-space crystallography concerns necessary ex- 
tinctions in the diffraction patterns of materials with certain non-symmorphic space groups. 
The simplest example of the n-point correlation function in (1), p(k)p(— k) = |p(k)p, is 
proportional to the intensity in a diffraction experiment. As is well known, materials with 
certain space groups show zero intensity at particular lattice points k. From the Fourier- 
space point of view, if some point-group element g leaves k invariant, and if ^g(k) 7^ 0, (3) 
requires that p(k) itself should vanish. 
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By (5), the condition gk = k means that all specimens $ in a particular equivalence 
class of phase functions will take the same value for that choice of g and k, regardless of 
gauge. Wc call the particular value ^^(k) a gauge-invciriant part. 

By definition, if two phase functions are related by a gauge (5), then they agree on all 
their gauge-invariant parts, but the converse is not true: out of the 230 periodic space groups 
in three dimensions, there are precisely two non-symmorphic space groups with no non-zero 
gauge-invariant $g(k). Thus they agree with their symmorphic cousins whenever gk = k, yet 
they are not symmorphic. Lacking non-zero gauge-invariants $g(k), the two space groups, 
/2i2i2i and /2i3, exhibit no systematic extinctions. For the other 228 periodic space groups, 
and for all known quasiperiodic space groups, it happens to be the case that two phase 
functions are related by a gauge if and only if they agree on all their gauge-invariant parts 
of the form ^^(k). If we knew this true in advance, we could use it to simplify classification 
of equivalence classes of phase functions by just finding all gauge-invariant parts compatible 
with (4). Absent this proposition, the gauge-invariant parts are insufficient: we must find a 
clever choice of gauge in which all non-invariant phases vanish. 

For the two peculiar cases in which one cannot make all non-invariant phases vanish, 
Mermin has constructed a gauge-invariant linear combination of two phase functions. 



for a specific choice of commuting point-group elements g, h and lattice vectors k/j and 
kg, where neither $g(k/j) nor ^fi(kg) is gauge- invariant. By showing that the group- 
compatibility condition (4) permits the gauge invariant (6) to assume a value either (sym- 
morphic) or 1/2 (non-symmorphic), he derives the existence of the two non-symmorphic 
space groups, /2i2i2i for point group 222 on the orthorhombic / lattice and /2i3 for point 
group 23 on the cubic / lattice. 

The fact that a gauge invariant of the form $g(k) is simply a linear combination of 
one phase, while (6) is a linear combination of two, suggests a generalization of the untrue 
converse about gauge- invariant single phases: 

Theorem 1. Two phase functions arc related by a gauge if and only if they 
agree on all their gauge-invariant linear combinations, of the form $g.(kj). 

We will show (Section 5) that an appropriate homology group systematically classifies all 
gauge invariants and (Section 6) that its duality to a cohomology group proves the theorem. 



Although a non-zero value of the invariant combination (6) implies no necessary ex- 
tinctions in the diffraction pattern, Konig and Mermin ■'^'^ have given it a definite physical 
interpretation. They consider a reciprocal-space vector q, not necessarily in C, and the 
associated little group, Gq C G, consisting of elements g such that 




(6) 



3. Necesseiry degeneracies 




(7) 
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is in the lattice. The Hamiltonian hq of the problem at Bloch wave- vector q describes all 
bands at q; bands may repel or cross. If for commuting g and h in Gq, the invariant (6) 
does not vanish, every level must be degenerate with some other level at q, so bands cross 
at least in pairs. We outline the argument, first because the gauge invariant (6) falls out of 
the first-homology calculation below and second because it relates to a second cohomology 
group. 

Konig and Mermin expand the Bloch wave in reciprocal-lattice vectors k and express 
the action of g in the point group on k in terms of the unitary operator U{g), with 

C/(^)|k) = e2-^«^'^)bk-k,) . (8) 

The operators U (g) commute with the Hamiltonian hq and constitute a ray representation 
of Gq, for 

U{g)U{h) = e-2'^^^«('"^)c/(5/i) . (9) 

The map 

\{g, h) = e-^'^^^^^^") (10) 

from G X G to the unit circle is called the factor system at q. 

Suppose that E is a non-degenerate energy level of hq, and let {ip) be the correspond- 
ing eigenvector. Since the Hamiltonian commutes with all the U{g), lip) must also be an 
eigenvector of each, so 

U{g)U{hm=U{h)U{gm . (11) 

Comparing this commutation of matrices U to the ray representation (9), Konig and Mermin 
find, for commuting g and h, that (6) must vanish. If it does not, every band must cross at 
least one other band at q. 

4. Gauge-equivalence classes and the first cohomology group 

Brown traces the theory of group cohomology to a 1904 work by Schur^^ on ray 
representations such as (9). We now use this theory to describe the classification of phase 
functions up to gauge equivalence. In Section 9, we compare this Fourier-space cohomology 
to group cohomology in real-space crystallography. 

The phase function ^y, which maps the lattice C to the real numbers modulo unity, 
has three defining characteristics: (i) it is linear, (ii) it satisfies the group-compatibility 
condition (4), and (iii) it is defined only modulo gauge functions of the form (5). Cohomology 
conveniently packages functions of just such description. 

The map $, which takes G G to is an example of a 1-cochain, the set of all 
of which we call C^. We denote by C* = Hom(£, R/Z) the group of all real-valued (mod 
1) hnear functions on the lattice £; the phase function is in jC*. ^'^ More generally, an 
n-cochain in G" is a function that maps n group elements from G to £*. has a group 
structure under addition. The gauge x, taking no point-group arguments, is a 0-cochain. 
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Connecting these functions is a coboundary map, d : ^ C""'"-'^, such that two maps 
in a row yield zero identically, d^'^'^^^ o d^"^^ — 0. Apphed specifically to the lowest-order 
cases, -"-^ 



(12) 



and {g, h) = ^goh-^gh + 



(13) 



Equation (12) identifies those 1-cochains that are coboundaries of 0-cochains as the 

d 

group of gauge transformations (by (5)); we define B'^ — im(C"' ~^ C") as the set of 
all n-coboundairies. Similarly, the right-hand side of (13) must vanish if $ is to satisfy 

the group-compatibility relation (4); we define = ker(C" ^ C""*" ) as the set of all 
n-cocycles, that is, n-cochains with vanishing coboundaries. 

The quotient group — jB^ ^ or cocycles modulo coboundaries, is called the first 
CO homology group of G with coefficients in L* . Membership in establishes group- 
compatibility, while modding out by S-*^ identifies functions differing only by gauge transfor- 
mation, so B^ is exactly the group of equivalence classes of phase functions. 



To prove Theorem 1, we must find all gauge-invariant linear combinations of phase 



functions to a real number; we may think of any gauge invariant in these terms. Since the 
value of a gauge invariant is the same for any member $ of an equivalence class we 
may write /({$}) = /($) as only a slight abuse of notation. We may turn this around to 
think of $ mapping / to a number. 

To produce a number in the unit interval, a cochain acts on integral linear combinations 
of ordered sets each containing one lattice vector and n group elements, which we can write 
(following Brown's convention -'^^) ^^i[fi'l«lfi'2il • • • Is'm] for kj G £ and gj^ E G. Since C 
already absorbs any integral coefficients, wc need not write them explicitly Such linear 
combinations are called n-chains and the set of all n-chains Cn = C[G'^'']. 

A cochain $ in G^ and a chain c = ^i[gi\ in Ci act on each other through the bracket 



We let Ci inherit the additive group structure of JC: 'k-lgi] + ^'[gi] = (k -|- k')[g'i]. 

By direct analogy with the coboundary operator in (12) and (13), we define a boundeiry 
map, 9 : Cn — > Gn-i, such that two maps in a row yield zero. Note that d decrements the 
number of copies of G in the chain, whereas the d operator on cochains increments it. In 
the lowest-order cases, 



5. Invariants and the first homology group 



functions. Suppose the linear function / : B 



R/Z takes an equivalence class of phase 




(14) 
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dHg] ^ gk-k 
and dk[g\h] = (hk) [g] - k[gh] + k[h] 



(15) 
(16) 



1-cycles are chains c = X^jkj[g(j] satisfying 



dc = ^{giki - ki) = 



(17) 



d 



Let Zi = ker(Ci ~^ Cq) denote the set of all 1-cycles. Applying the bracket (14) to c 
satisfying (17), we find that 1-cycles give gauge-invariant linear combinations of phases, and 
since the condition for gauge invariance of a general linear combination of phases <l>g.(kj) 
is, by (5), precisely (17), all such gauge-invariant linear combinations come from 1-cycles. In 
the case that c = k[g] in (17), the gauge invariant takes the form ^g{k) for gk — k, but in 
general a linear combination, such as (6), must be allowed. 

Similarly, l-boundeiries are those 1-chains that can be written as boundaries of 2-chains, 

d 

as in (16) or integral linear combinations thereof; Bi — im(C2 ~^ Ci). Because d o d = 0, 
every 1-boundary is a 1-cycle, representing a gauge invariant. Again applying the bracket we 
discover that the phases at boundaries (16) are necessarily zero by the group-compatibility 
condition (4). Thus, although 1-boundaries are gauge invariants, they are trivial. We mod 
them out in defining the first homology group of G with coefficients in C: 



First homology is the group of non-trivial gauge invariants. 

We have already commented that a real- valued function f on is unaffected by gauge 
{i.e., by the addition of a coboundary to a cocyle $) and that the addition of a boundary 

to a cycle x does not change any phase. These facts establish that the bracket, which we 
introduced in (14) function x Ci ^ R/Z, is in fact well defined on x Hi, which 
we may express concisely thus: 



where $ is a phase function (rf$ = 0), a gauge transformation of the form (5), a; G Ci a 
cycle {dx = 0), and dy a boundary, y e C2- 



We note a striking similarity between H and Hi. While the former contains all linear 
functions satisfying group compatibility, identifying those related by gauge transformation, 
the latter contains all linear gauge invariants, identifying any whose difference is made trivial 



Hi = Zi/Bi 



(18) 



($ + dx,x + dy) = ($, x) 



(19) 



6. Duality and the Proof of Theorem 1 
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by group compatibility. In fact, homological algebra affirms their duality. Before stating this 
central result, we note the following: 

Lemma 2. For G a finite group with \G\ elements and $ e H^, \G\^ = (where 
G denotes the set of phase functions gauge-equivalent to zero). In other words, 
for any phase function $, there exists a gauge in which every $g(k) is rational with 
denominator \G\. 

This is quite plausible, since the only constraints on phases are those imposed by group 
compatibility (4) applied through the constituent relations of the point group. Specifically, 
summing over g in (4) and noting that ^gf^ = we find 

|G|$;,(k)=x(^k-k) , (20) 

where x(k) = — Ylg ^gO^) clearly satisfies the definition of a gauge function. Then |G| times 
any phase function is a pure gauge, proving the lemma. Note that we may not divide by 
|G| to show that all phase functions are gauge-equivalent to zero, because x is defined only 
modulo unity. 

Having found a gauge in which a phase function takes rational values, we define C' — 
Hom(>C, (Q/Z) C jC* and quote the following standard result from cohomology: 

Lemma 3. For G a finite group and C a free Abelian group, the groups 
//"(G, £') and Hn{G, C) are dual through the bracket (14) (,)://" x ^ Q/Z. 

Duahty means that if X serves as a minimal basis for ifi, there exists a basis {$^} for 
labeled hy x & X and that any $ G can be decomposed uniquely as a linear combination 
of the basis elements <I>^ with integral coefficients equal to ($, a;)/($^, x). Therefore, if 
we know the value of a phase function on each of the non-trivial gauge invariants x 

spanning Hi^ by linearity we know $ up to gauge transformation. This proves Theorem 1. 

We give in the appendix an alternative, elementary proof of Theorem 1 by constructing 
the gauge function relating two phase functions that agree on all gauge-invariant linear 
combinations. 

7. Calculating Space Groups (example) 

Homological algebra offers a powerful tool for systematically finding all gauge invariants. 

By Theorem 1, it furthermore frees us from having to specify a gauge: once we have the 
gauge invariants, we no longer have to prove on a case-by-case basis that all other phases 
can be made to vanish. As an example, we compute the two possible space-group classes 
on the body-centered orthorhombic direct-space lattice (which is face-centered in reciprocal 
space) with the point group 222. Instead of having to hunt for the gauge invariant (6), we 
shall see it fall out of the formahsm naturally. 

The calculation is straightforward. First we calculate the group Zi of gauge invariants by 
acting with the boundary operator d on 1-chains \^[g\. We then wish to remove from the list 
all boundaries in Si, since these represent the "trivial" consequence of group compatibility 
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(4); to do so, we let d act on 2-chains kfo^l/i]. After this elimination, we find a homology 
group with only two elements. The entire procedure is easily automated; for example, the 
NullSpace[] and LatticeReduce [] functions of Wolfram's Mathematica program can be 
used to find minimal bases over Z for Zi and Bi. 

Let e*, i = 1,2,3, constitute the usual Cartesian axes. We name the three sides of the 
conventional reciprocal cell x = ae^, y = be^, and z = ce^, where a, b, and c are all different, 
and generate the face-centered reciprocal lattice C (dual to the body-centered direct lattice) 
with the vectors bi — |(011), 62 = ^(lOl), and 63 = ^(110) (all in terms of the (x, y, z) 
basis). The point group is G = {e,di,d2,ds}, where dj leaves e^' invariant, did2 — d^ et 
cycL, and df — e, the identity. G is generated by di and ^2- 

To calculate the group of gauge invariants, Zi, we apply d to the (additive) generators 
of 1-chains. The chains k[e] are cycles, since e leaves all vectors invariant. We can show that 
they are also boundaries, for by group compatibility, 

k[gh] = {hk)[g] + k[h] , (21) 

where = denotes equality up to the boundary ^kfgfl/i]. Setting g ~ h — e, we then have 
k[e] = 0, so we leave it out of the computation. Equation (21) furthermore tells us that we 
need consider only the generating elements di and d2, not ^3 — did2- 

One easily verifies that 

{dhi[dj])k = \eijk\ - - Sjk , (22) 

where the subscript k indicates the resulting lattice component in the (bi, b2, h^) basis and 
where 6 is the Kronecker delta and e the totally antisymmetric tensor. 

With i = 1,2,3 and j restricted to 1,2, we can think of the left-hand side of (22) as a 
thrcc-by-six matrix whose three-dimensional null space over the integers is Zi (aside from 
the boundaries we have already eliminated). It is a standard result of linear algebra that 
row reduction without division yields a primitive basis over the integers; we find for this 
basis 

zi = x[di] 

Z2 = y[d2] (23) 
Z3 = -bi[d2]+b2[di] . 

We now know all the gauge invariants and seek, by computing Bi , to identify those that are 
related by group compatibility (4). 

To find the 1-boundaries not already identified, we must apply d to the 27 two-chains 

hi[dj\dk] ; (24) 

this number is reduced to 18 if we notice that k[g'|/ii/i2] has the same boundary as (/i2k)[5'|/ii] 
-I- k[g'/ii|/i2] — k[/ii|/i2], so that in (24) we may restrict A; to 1,2. 
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After eliminating duplications and zeroes, we get seven boundaries. Rearranging the 
rows into echelon form, we find only three that are linearly independent: zi, Z2, and 2z^ 
are all boundaries. This identifies all even multiples of the cycle zs in Hi, which therefore 
contains just a single non-zero invariant: 

HiiG,jC)^{0,Z3} . (25) 

Since 2z^ = 0, $(-23) can take values only and 1/2. 

The invariant is precisely the linear combination (6) of phases used in Ref. [5] to 
distinguish the symmorphic (/222, $(-23) = 0) from the non-symmorphic (/2i2i2i, $(,23) = 
1/2) space group. Mermin's rather shorter calculation rests on a clever choice of gauge 
ensuring ^^^(h^) = ^^^(hi) = ^^^{102) = 0; in our comparatively pedestrian approach, 
the list of phases that can be discarded inheres in the calculation of the boundaries Bi. 
Significantly, the present calculation is easy to automate, requiring no "judicious choice" of 
generators and gauges, while verification that one gets "a representative of every class" ^ is 
built in to the procedure. We lose the elegance of Mermin's (equivalent) calculation, but we 
hope we make up for some of the loss in the compact and general statement of Theorem 1 
as proven in Section 6. We have generalized and applied essentially the ideas demonstrated 
in this section to the simultaneous computation of a large number of cases, including non- 
standard and modulated lattices. 

8. Second Cohomology and the Factor System 

In general, the cohomology group H'"'{<3,dJl) consists of cocyclcs from n copies of a 
group C5 to a module 971 with the identification of coboundaries. Having classified the phase 
functions in with <3 = G and OJt = C, wc can calculate the factor system (10). However, 
the factor system also has a cohomological existence in its own right, with a different group 
<3 and different module 971. Consider 

A(^,/i) = -(logA)/(27ri) = $g(k,,) , (26) 

taking values in the module 971 = [0, 1) = R/Z. From here on, we will use the term "factor 
system" to refer to A rather than to A. Since g and h live in the httle group = Gq, A is a 
2-cochain in C^(Gq, R/Z). The action oi g E Gq on x in the module R/Z is trivial: xg — x. 

We can again impose group compatibility as a cocycle condition: 

{d^){9l, 92, 93) = H91, 92) - H9l,929s) + H9192, 9s) - H92, 93) ^ ^ ■ (27) 

(The equivalence to group compatibility follows from the definitions (7) and (26).)^^ 

If two ray representations (8) differ by a quantum phase e^""^^"^*^^) (one for each g e Gq) 
so that U'{g) = e^'^'^'^^3^U{g), the corresponding factor systems differ by the coboundary of 
(t: 

K{g,h)-h'{g,h)^{5a){g,h) = a{g)-a{gh) + a{h) . (28) 
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Since the quantum phase has no physical consequence, we identify ray representations re- 
lated by (28). The set of phase- equivalence classes of ray representations is the second 
cohomology group, H'^{Gc^, R/Z) = Z'^/B'^. 

As a final application of group cohomology, we consider a lattice vector k left invariant by 
a point-group element g & G. Letting & = (g) = {g"'\n G Z}, the cyclic group generated by g, 
and 971 = £, wc express k G £ as a 0-cochain. The cocyclc condition {d'k){g"') = g^k — k = 
expresses the requirement that k be fixed by g. Since there are no (— l)-cochains, there is 
no identification of coboundaries, and the cohomology — is the group of all k left 
invariant by g. 

9. Comparison to Real-Space Cohomology 

It has long been known that cohomology groups can be used to classify the space groups 
of periodic crystals; Hiller^^ gives a particularly approachable exposition. We now compare 
our formulation, using the Fourier-space lattice JC, to the traditional approach, described in 
terms of the periodic lattice L in real space. In terms of L, we can describe jC as the lattice 
of integer- valued homomorphisms on L, C — Hom(L, Z). (For mathematical convenience, 
and to allow direct comparison with Hillcr, we have in this section absorbed the conventional 
factor of 271 into C.) In the three-dimensional case, if we choose a basis ai, a2, as of L then 
we can describe C in terms of the dual basis ki, k2, ks, where (kj|aj) = S^j (Kronecker 
delta). 

Any real-valued linear form on C has the form (k|x) for a fixed vector x in real space. 
Recall that our phase functions $^ take values in R/Z, not R: that is, we identify ^^(k) 
with $g(k) + n, for any integer n. This corresponds to identifying x with x -|- a, for any 
a G L. Thus there is a one-to-one correspondence between phase functions and elements of 
the quotient group R^/L, and we can write $g(k) = (k|xg) for some real-space vector x^, 
identified with any of its translates by lattice vectors. 

In Fourier-space crystallography, the phase functions are required to satisfy the group- 
compatibility condition (4), or equivalently the cocycle condition {d^){g,h) = (c/. (13)). 
In terms of the vectors x^ in real space, this condition becomes (k|xg/j) = (/ik|xg) + (k|x/j). 
Recall what it means to let h (an orthogonal transformation on real space) act on k (a vector 
in Fourier space): by definition, (/ik|x) = (k|/i~^x). Thus the group-compatibility condition 
becomes 

^gh = h~\ + ■ (29) 

In order to compare this with the notation of Ref. [24], let s{g) = x^-i and replace g and h 
with g~^ and h~^, respectively: 

s{hg) = x^-i^-i = x/j-i + hxg-i = s{h) + h s{g) . (30) 

This is exactly the cocycle condition used in Ref. [24], Prop. 5.1. 

Finally, consider gauge equivalence: we identify phase functions and related by a 
gauge transformation as in (5). In the periodic case we are considering, the gauge function 
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can be described by x(k) = (k|y), where again y is a real-space vector considered modulo L. 
If the second phase function $^ corresponds to in real space then the condition (5) can 
now be expressed as (k|x^) — (k|X(^) = ((/k — k|y), or x^ — = g~^y — Y- In terms of 
s{g) = x„-i and s\g) = x' this becomes 



which is exactly the coboundary condition in Ref . [24] . 

In conclusion, the cohomology group = Z^/B^ that we consider here agrees, in the 
periodic case, with the cohomology group H^{G, R^/L) described in Ref. [24]. 
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We have already shown that duality proves Theorem 1, which wc may restate, a phase 
function is gauge- equivalent to zero if and only if it vanishes on all gauge-invariant linear 
combinations, or again cquivalently, a 1-cocycle is a coboundary iff it vanishes on all cycles. 
The implication ( "only if" ) statements come immediately, so we prove the converse. 

Although standard treatises on homological algebra prove duality of and Hn (Lemma 
3) , the following elementary demonstration of Theorem 1 also explicitly constructs the gauge 
transformation x for which $^(k) = x{9^ — k) if $ vanishes on all gauge invariants. The 
proof extends familiar ideas from the linear algebra of vector spaces to that of free modules 
(over Z). 

Define the linear function x on all 0-boundaries do E Bq {c E Ci) by 



We can do so consistent with linearity because, by hypothesis, ($, z) = for all gauge 
invariants z e keid. As a special case of (32), for c = \^[g], we have ^^(k) = x{9^ ~ k), 
which almost makes $ gauge-equivalent to zero. However, to establish that $ is purely a 
gauge (thus proving the theorem), we must extend x, which we have so far defined only 
on Bq, to a linear function x defined on all C — Cq. To extend x to x on Cq, write the 
rectangular matrix whose Uc = rank(So) columns are the integral expansions of a minimal 
basis of Bq in a basis of Cq; M has rir — rank(Co) ^ rows. M maps Bq to Cq by inclusion. 

M M is a Uc X Uc matrix; if M M is not singular, we can form the left inverse of 
M as M = {M^M)-^M^. To show M^M non-singular, observe that the columns of 
M are all independent over the integers, so that we can find a matrix M' with Uj- rows 
and n,. — columns, all columns orthogonal to those of M . The x Ur square matrix 
A = [MM') is therefore non-singular. Since A and A^ are both non-singular, so is A^A = 



s{g) - s'{g) ^y-gy 



(31) 
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Appendix: Elementary Proof of Theorem 1 




(32) 




, but since det(^^ A) ^ 0, we must have det(M^ M) ^ 0. 
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Think of x as a row vector of rank nc] define x — (row vector of rank rir)- This 
completes the construction of a hnear x on all Cq such that ^^(k) = x(5'k — k), proving 
Theorem 1. 
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